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We investigate the late-time behavior of the massive vector field in the background of the 
Schwarzschild and Schwarzschild-de Sitter black holes. For Schwarzschild black hole, at interme- 
diately late times the massive vector field is represented by three functions with different decay 
law * ~ t~ {t+ ' i/2) sinmt, *i ~ t~ < - l+5/2) smmt, * 2 ~ i~ (<!+:L/2) sin mi, while at asymptotically late 
times the decay law $ ~ t~ 5 ^ 6 sin (mt) is universal, and does not depend on the multipole number I. 
Together with previous study of massive scalar and Dirac fields where the same asymptotically late- 
time decay law was found, it means, that the asymptotically late-time decay law ~ t~ 5 ^ 6 sin (mt) 
does not depend also on the spin of the field under consideration. For Schwarzschild-de Sitter black 
holes it is observed two different regimes in the late-time decay of perturbations: non-oscillatory 
exponential damping for small values of m and oscillatory quasinormal mode decay for high enough 
m. Numerical and analytical results are found for these quasinormal frequencies. 
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I. INTRODUCTION 

Black hole's response to an external perturbation has 
been a subject of active investigations for recent ten years 
[l[ . Now it is a well known fact that in the response sig- 
nal of a black hole, the so-called oscillatory quasinormal 
frequencies dominate, which do not depend on excita- 
tions of the quasinormal modes, but only on a black hole 
parameters, thereby giving us the complete information 
about the geometry of a black hole in the fully non- linear 
general relativity theory [2[ . 

As was shown for the first time by R. Price, the regime 
of quasinormal oscillations turns into decay with inverse 
power dependence on time for a Schwarzschild black hole 
[3[. Price showed that perturbations of the massless 
scalar and gravitational fields decay as i~( 2f + 3 ) at asymp- 
totically late times t — > oo. Bicak found that scalar mass- 
less field in the Reissner-Nordstrom background decays 
as i-( 2£ + 2 ) for \Q\ < M, and as i~^+ 2 ) for \Q\ =M|. 
The massless scalar field perturbations along null infinity 
and along future event horizon has the decay law u~ l+2 
and v~ e+3 respectively, where u and v are outgoingand 
ingoing Eddington-Filkenstein coordinates Q. In Q it 
was shown that the charged scalar perturbations decay 
slower then a neutral ones at asymptotically late times, 
while, on the contrary, at the stage of quasinormal ring- 
ing, the neutral perturbations decay slower @, §• Also, 
for Schwarzschild-de Sitter and Reissner-Nordstrom-de 
Sitter black holes, instead of power-law tails, the expo- 
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nential tails were found 0]. 

The investigation of late time decay for massless fields 
was continued for Z?-dimensional black holes. Thus, 
tails of massless scalar, vector and gravitational fields 
for higher-dimensional Schwarzschild black holes have 
the decay law £-( 2 ^+-°~ 2 ) for odd D > 3, and for even 
(D > 4) dimensions, the power law decay is £-( 2£ + 3 -°- 8 ) 
The late-time behavior for £>-dimensional Gauss- 
Bonnet, Gauss-Bonnet-de Sitter and Gauss-Bonnet-anti- 
de Sitter black holes was elaborated in Late time 

tails were also observed perturbations around wormholes 
between the branes in a two brane Randall-Sundrum 
model [Hj ■ In addition, the late-time tails can be a tool 
to probe the extra dimensions both at the ringing stage 
[l3l | and at the stage of late-time tails • 

The late time behavior of massive fields is qualitatively 
different from massless ones: at late times the decay pro- 
file is oscillatory inverse power tail. For a massive scalar 
field with mass m in the background of the Schwarzschild 
black hole with mass M, the perturbation decays as 
£-(^+3/2) s i nm f at intermediate late times mM < mt < 
1/ (m M) 2 I16l . and as £~ 5 / 6 sin mi at asymptotically late 
times The same ~ i -5 / 6 behavior was found for 

the massive scalar field perturbations of the Kerr black 
hole [l9[ (but with frequency which is slowly dependent 
on time ^ ~ i~ 5 / 6 sin (u>(t) x i)), of the dilaton black 
hole [2(J, and also, for massive Dirac perturbations of 
the Schwarzschild black hole [2l[. Recently, it has been 
found in [22| that intermediate late time decay of massive 
scalar field in D-dimensional Schwarzschild black hole is 
proportional to t e+D / 2 ^ 1 / 2 . 

On the other hand, the late time tails of fields corre- 
sponding to massive bosons have not been studied before, 
even for the simplest and a most interesting background 
geometry, namely, for the Schwarzschild black hole. In 
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the present paper we shall investigate the late time be- 
havior of massive vector field, which shows quite non- 
trivial behavior also at the stage of quasinormal ringing 
f24l ]. The massive vector field obeys the Proca equation 
which can be reduced to a single wave-like equation for 
spherically symmetrical perturbations [24[ . Therefore we 
are able to investigate the tail behavior for spherically 
symmetric perturbations both numerically and analyti- 
cally. Yet, for higher multipoles £ > the Proca equa- 
tions cannot be decoupled. Fortunately, in this case we 
also can investigate the tail behavior, because we need 
only the asymptotic form of the wave equations at large 
r, rather than its exact form. 

We found that at intermediately late times the three 
types of perturbations of the vector field decay with three 
different decay laws, while, at infinitely late times, these 
three types of decay approach a single universal decay 
law, which is independent of multipole number £, and is 
the same for massive fields of other spin. Late time be- 
havior for asymptotically de Sitter geometries essentially 
depends on the value of the mass of the field m: at some 
treshold value of m the non-oscillatory decay changes into 
the oscillatory one. 

The paper is organized as follows: in Sec. I we study, 
both analytically and numerically, the late time behavior 
for spherically symmetrical perturbations of Proca field 
in the background of Schwarzschild black holes. Sec II. 
is devoted to the analytical consideration of late time be- 
havior for non-spherical perturbations. Finally, in Con- 
clusion we summarize the obtained results. 



II. SPHERICALLY SYMMETRIC 
PERTURBATIONS 

We shall consider here the Schwarzschild black hole 
solution with a positive cosmological constant A, i.e. 
Schwarzschild and Schwarzschild-de Sitter backgrounds 
in which the massive vector field propagates. The black 
hole metric is given by 



where 



-h(r)dt 2 + hir)- 1 ^ 2 + r 2 (d0 2 + sir 



h(r) = 1 - 



2M 



Ar 2 
~3~' 



; ), (1) 
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For self-consistence, let us briefly deduce the wave 
equation for spherically symmetric perturbations of the 
massive vector field for general form of the spherically 
symmetric static metric (1). The Proca equations have 
the form: 

Fg - m 2 A» = 0, F^ = A Utfl - A^ v . (3) 
These equation can be reduced in the spherically sym- 



metric case to the single equation of the form: 
B,tt , (2f(r) 



f(r) V r 



+ f'(r) B, r + 



/2/'(r) 2/(r) 
\ r r 2 



m 2 ) B = 0, 



where 



B = A 



r,t At t r 



(4) 
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Introducing the usual "tortoise" radial coordinate r + 
and a new wave function 

*(r) = B(r)r , (6) 

the wave equation can be re- written in the form [2^ |: 

d 2 ^{r+,t) d 2 V(r+,t) 



dr 2 dt 2 
with the effective potential 



V(r(r*))*(r*,i) = 0, (7) 



V(r) 



1 



2M Ar 2 



6M 



(8) 



Positive values of the cosmological term A correspond to 
asymptotically de Sitter solutions. 

The time evolution of the massive vector field can be 
described by the spectral decomposition method [23|. 
Thus, 

*(r*,t) = 

= y'ri(G(r*y i ;t)*tW,0) + G t (r*, 0)). (9) 

For t > 0, the Green function obeys the equation 
d 2 G(r*,<;t) d 2 G(r*,<;t) 



dr 2 



8t 2 

6(t)8(r* O 



V{r)G{r*y^t) 



(10) 



In order to find the Fourier's component of the trans- 
formed Green function G (rv,r^;t), we need the two lin- 
early independent solutions '5 2 to the homogeneous 
equation 



dr 2 



(u 2 -V)Vi = 0, i = 1,2. 



(11) 



If A = 0, we require that the function vPi should rep- 
resent purely in-going waves at the event horizon (7' = 
r + ): ^1 ~ e~ twr * as r ^ _ > —00. The other func- 
tion \P 2 should damp exponentially at spatial infinity: 
v|/ 2 r~j e -Vu' 2 -m 2 r* ^ ag r ^ _^ ^qq^ Q n other hand, if 
A > 0, the spacetime have a new feature, the cosmolog- 
ical horizon (r = r c ). In this case, the region of interest 
is the block r + < r < r c . The effective potential decays 
exponentially to zero at both limits, and therefore we 
require that ^2 ~ e~ luJT * , as r* — + +00. 
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To check the analytical results which will be found 
below, we used a numerical characteristic integration 
scheme, based in the light-cone variables u = t — and 
v = t + T*. In the characteristic initial value problem, 
initial data are specified on the two null surfaces u = uq 
and v = vo- The discretization scheme applied, used for 
example in [25[ , is 



0(A 4 ) 



(12) 



where we have used the definitions for the points: N — 
(u + A, v + A), VF = (it + A, u), £ = [u,v + A) and 
5 = (it, u). 

In the asymptotically de Sitter scenario, we also em- 
ploy a semi-analytic approach to investigate the pertur- 
bative behavior of the vector field. As will be seen, there 
is a region of the parameter space of the geometry where 
the late-time perturbations are dominated by quasinor- 
mal mode behavior. At this regime, we use the WKB 
formula 



Vo 



hi — La 



(13) 



where Vq is the height of the effective potential, and Vq 
is the second derivative with respect to the tortoise coor- 
dinate of the potential at the maximum. L2, L3 L4, L5 
and Lq are presented in [l5|. 



A. Schwarzschild black hole background 

1. Intermediate late time behavior 




Figure 1: Late-time behavior of the massive vector field per- 
turbation (Mm = 1) in asymptotically flat geometry, at 
r = 4.5M (initial conditions have been chosen constant = 1 
on the null cone). The analytically obtained law ^ oc t~ 5 ^ 6 
is shown by the line. At late times numerically we have 
the power coefficient —0.844, what differs from —5/6 only 
by 1.3%. 




At intermediately late times, when r <C t <C 
M/ (Mm) 2 , one can neglect the effect of backscattering 
of the field from asymptotically far regions and, following 
[ill, use the approximation M Cr M / (Mm) 2 . One 
can expand the wave equation into power series in M/r. 
Neglecting terms of order 0((Mm/r) 2 ), we obtain 



dr 2 



AMuj 2 - 2Mm 2 



4) ! (1-4! 



Provided that the observer and the initial data are situ- 
ated far from the black hole we can further approximate 



d 2 ^ ( 2 2 2 
u> — m 



dr 2 



* = 0. 



(15) 



The last equation is similar to Eq. (15) in [16J, therefore 
we can easily find the branch cut contribution to the 
Green function for the above equation: 



9ttA 2 



G C (r„r:,t) = 
yb((r+,k)yb((r'„k)(k) 3 e- 



'dio 



(16) 



Figure 2: Late-time behavior of the massive vector field per- 
turbation (Mm = 1.5) in asymptotically flat geometry, at 
r = 1.6M (initial conditions have been chosen constant $ = 1 
on the null cone). The analytically obtained law ^ oc i~ 5//6 
is shown by the line. At late times numerically we have 
the power coefficient —0.849, what differs from —5/6 only 
by 1.76%. 



where k = \/m 2 — lo 2 , and A is a normalization constant. 

Note that the contribution to the Green function is 
dominated by low-frequencies at late times. In the large 
t limit the effective contribution to the above integral 
is from y/m 2 — uj 2 = O(ymjt). This is stipulated by 
rapidly oscillating term e~ iujt so, that there is cancella- 
tion between the positive and negative parts of the inte- 
grand. Finally, similar to [161 ], for t 3> m^ 1 we obtain for 
a fixed radius, and r*, <C t: 



2 m 3 / 2 (r^<) 2 t- 5 / 2 cos (mt - (5?r/4)) 



6! 



G C (r*,rl,i) = 

V ,, u: 

Therefore one can conclude that at intermediate times 
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the field \& decays at a fixed radius according to the law 
* ~ t~ 5/2 smmt, (18) 

and at the black hole event horizon, the decay is the same 
$ ~ i,- 5 / 2 sin mt. (19) 

This analytically obtained result in the regime of inter- 
mediately late times is confirmed by numerical character- 
istic integration with high accuracy, what is illustrated in 
Figs. [T]and [2 There one can see that at sufficiently late 
times the numerical envelope approaches the analytical 
law \I/ oc i -5 / 6 . The difference between numerical enve- 
lope and the line t -5 / 6 is decreasing as time increases and 
is about 1% for t ~ 10000 (see Fig. QJ. 

2. Asymptotic late-time behavior 

In (l7j , it was shown that at asymptotically late times 
the law of decay is stipulated by the behavior of the wave 
equation at large r. Thus, in the region r/M 3> 1, ex- 
panding the wave equation in powers of M/r, we obtain 



dx 2 



1/4 



V = o, 



(20) 



where we used a new variable x = 2r\Jm 2 — uj 2 , and a, 
b are constants depending on m, M, and lu. Following 
the papers [l8[ and [13] , one can show that the monopole 
massive vector field at asymptotically late times t — > oo 
undergoes oscillatory inverse power law decay 



* ~ t 5/6 sin mt, t -> oo. 



(21) 



This power-law envelope in shown in Figs. [T] and [5J In 
the range mt l/(mM) 2 , the smaller the value mM is, 
the later the i~ 5 / 6 tail begins to dominate. In the range 
mt 3> mM, on the contrary, the larger the value Mm is, 
the later t~ 5 ^ 6 tail begins to dominate. The numerical 
investigation show that, both in the intermediate and 
late time regime, the oscillation period can be perfectly 
approximated by 



period 

This is illustrated in Fig. [3[ 



2tt 



(22) 



B. Schwarzschild-de Sitter black hole background 

The introduction of a positive cosmological constant 
changes drastically the late-time tails. Direct numerical 
integration shows two distinct late-time decay regimes. 
We have observed that the qualitative aspects of the field 
decay are independent of the value of the cosmological 
constant, in the non-extreme regime. 
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Figure 3: Numerical value of t in which the wave function is 
zero. Strait lines imply that period of oscillations is a con- 
stant. For late times, the lines t zero (i) = 2ni/m are excellent 
approximations of the numerical values. 
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Figure 4: Semi-log graphs of the massive vector field pertur- 
bation in de Sitter geometries, at r* = 0. Increasing the mass, 
it is observed the transition from non-oscillatory to oscillatory 
asymptotic decay. In the graphs, M = 1 and A = 10 -4 . 



For fields with small values of mass, we observe that 
the late-time decay is exponential and non-oscillatory: 



(23) 



The exponential index k\ approaches zero as the field 
mass tends to zero. This decay mode is similar to 
the tail observed in the massless vector perturbation in 
Schwarzschild-de Sitter black hole backgrounds. 

For m greater than a critical value, the decay is domi- 
nated by quasinormal modes. In this regime, it is possible 
to use the well known semi-analytical WKB approach to 
calculate the quasinormal frequencies. In this approach, 
by the way, the convergence is very good. We present 
some results in table IIIB[ where it is compared frequen- 
cies calculated directly from the characteristic integration 
and first and sixth order WKB formulas. 

The transition from non-oscillatory to oscillatory decay 
is shown in Fig. [4] We also observe that, the larger the 
cosmological constant, the greater is the critical mass m 
where the oscillatory regime dominates. For high values 
of to, the exponential envelope index ki approaches a 
constant value. These observations are illustrated in Fig. 

El 



5 



Table I: Fundamental quasinormal frequencies for the massive vector field perturbation around the SdS black hole for several 
values of the parameters m and A. The black hole mass is set to M — 1. 

Direct Integration WKB-l" 1 order WKB-6" 1 order 
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Figure 5: Graphs of ki and ki with the mass m. The verti- 
cal dotted lines approximate the masses where it is observed 
the transition from non-oscillatory to oscillatory asymptotic 
decay. In the graphs M = 1. 
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Figure 6: Graph of Re(o;o) 2 as a function of m 2 , in the near- 
extreme limit positive A limit. The bullets are the values 
calculated from the time-evolution profiles, squares are values 
obtained from WKB formula, and the strait line represents 
the results from the expression (]24[) . The parameters in this 
graph are A = 0.11082 and M = 1. The differences between 
the results are under 1%. 



Following the work presented in [26|, |27( , we have 



One interesting limit of the Schwarzschild-de Sitter 
black hole geometry is its near extreme limit, when 
the event and cosmological horizons are very close, and 
the cosmological constant approach its maximum value 
(1/9M 2 ). 

In the near extreme limit, an analytic expression can 
be obtained for the quasinormal frequencies, which dom- 
inate the late-time decay if the mass m is large enough. 



Re(w n ) 



A(l - 9M 2 A) jm 



A 4 ' 



T . , A(l - 9M 2 A) / 1 
ImK) = !s '- f n+ - 

n = 0,l,2,... . 



(24) 



The above formula is very well confirmed numerically: 
Figl6]shows a comparison of the values obtained by direct 
numerical integration, WKB formula and from Eq. (I24p . 



III. PERTURBATIONS OF HIGHER 
MULTIPOLES 

In this section we shall consider only the case of 
Schwarzschild black hole geometry. It is known from the 
work of Galtsov, Pomerantseva and Chizhov [28[ that 
the perturbation equations for massive vector field in the 
Schwarzschild background cannot be decoupled for gen- 
eral case of arbitrary multipolarity I. Yet, far from black 
hole the equations can be decoupled |28| , what allows us 
judge about asymptotically late-time behavior. As in the 
paper [28] the perturbation equations are not shown in 
explicit form, but rather as an operator equations , w e 
shall have to repeat here some of the derivations of [2a ]. 
in order to proceed our analysis of tail behavior. 

The Proca equation reads 



F$ v -m 2 A^ = 0, F^ = A v , 



_ A 

[1 S^-fJ,,!/' 



(25) 



Using the Newman-Penrose formalism, we go over to the 
tetrad projections of the field 



A m . = J2 e-^ t+im *UY em )^{ m (r), (33) 
t=i,\m\<e. 



An= Yl e~ luJt+tm HoY em )y e 5 m (r), (34) 
Ai= Yl e-^ t+im ^( ^ m )^ 6 ro (r), (35) 

£=0,\m\<£ 

one can exclude the functions ^3, ^4, ^5, ^>q and obtain 
the three equations for the three scalar functions ^0, ^1, 
and * 2 HI: 

(D Q D+ - C/)* + \/2Ar- 2 *x = 0, (36) 



(D+D - C/)* 2 + v^Ar- 2 *! = 0, (37) 



*o = ^!V, $! = -F^ l n v + m^m v ), (26) 



= F, w m*»n\ 



(27) 



"2 — * JLII/' 

and of the vector-potential 

A" = Ain" + AJ" - A m rn*" - A m , m ". (28) 

After projecting of the equations (15) onto the tetrad's 
vectors and expansions of <f>i and A a into spin-weighted 
spherical harmonics [30l | 



r 3 (D+r- 4 D r + D 1 r- 4 D+r)V 1 - 2U^+ 

+2V2(r 2 -2Afr)- 1 (*o + * 2 ) = 0. (38) 



Here we used the operators Dj, j = 0, 1: 

iwr . 2(r — M) 
3 ~° r r-2M +J r(r-2M) 



(39) 



*o = e- iwt+im *d^ m )^ m (r) 

£=l,\m\<e 



(29) 



, „ icur 2(r-M) . , 



$1 = 



OO 

e=o,\m\<e 



- iut+im *(oY em )^{ m (r), (30) 



Here X = £(£+ 1), and 

U = (A 2 + ^V)(r 2 - 2Mr)~ 1 



$2 = e- iut+im ^UY im )^fi m (r), 
t=l,\m\<t 



A m = e-^'+^iY^^Cr), 
e=i,\ m \<e 



(31) 



Using new notations: 



x = I, u = 2V2mM, w = 2Mw 

2M r 



(32) we obtain the system of equations 
I 



(41) 
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x 2 x(x + 1) 
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(42) 
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*2,x* + 
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A 2 



x 2 x(x + 1) 



* 5 



V2A 
x(a; + 1) 



*i = 



(43) 



x(x + 1) 



+^^ TT (*o + * 2 ) = o. 

X(X + 1) 



A 2 



x- 2x 2 



x 2 (x + l) 2 x(x + l) 



(44) 



To check these equations one can consider here spher- 
ically symmetric perturbations, for which "Jo = ^ = 0, 
while v^i is non-zero and obeys the wave equation (|44p 
with I = 0. After coming back to r-coordinate, Eq. (|44|) 
reduce to Eq. ©. 



A. Intermediate late-time behavior 

Now, let us consider phcnomcnologically interesting 
region mM <C 1, what corresponds to the notion of 
"small" perturbations, i.e. perturbations where the mass 
of the field m is much smaller then the mass of the black 
hole. Again, at intermediately late times (r <C t <C 
M/(Mm) 2 ), we neglect the effect of backscattering of 
the field from asymptotically far region, and use the ap- 
proximation M <C r <C M l '(Mm) 2 . For an observer and 
an initial data situated far from the black hole, from the 
above Eqs. (|4"2")) - (|4"4")) . we obtain the following equations 



^0,xx 



*2 



t> 2 


A 2 > 




V2A 


2 


X 2 J 


x z 




A 2 ^ 


*2^ 


V2X 


2 


X 2 ) 


2 



*1, 



2 m 2 A 2 + 2 



*i = (45) 
Y&i = (46) 

/2A 

^V(*o + * 2 ) =0. 

The above equations can be easily diogonalized by lin- 
ear x-independent transformations of the ^-functions. 

Vo,xx +(w 2 ~^-^)yo = (48) 

* 2 = (49) 



2,xx 



EL 

2 



1 + A 2 - Vl + 4A 2 ^ 



^l,xx 



2 M 
u; 

2 



1 + A 2 + x/ITlA 2 



*i = (50) 



The above equations can be re- written in the following 
way: 

^O.xx+ W 2 



*o = 0, £ = 1,2,3,. 



(51) 



^>2,xx + 



2 M 2 ^ - 1) 



* 2 = 0, £=1,2,3,... 

(52) 



2 M 

ur 

2 



£(£ + 3) + 2 



*i=0, £ = 0,1,2,. 



(53) 

One can see that £(£ +1) in Eq. (|5Tj) runs val- 
ues 2,6,12,20,..., £(£ - 1) in Eq.© runs the val- 
ues 0, 2, 6, 12, 20, . . ., while £(£ + 3) + 2 runs values 
2, 6, 12, 20, .... Therefore we can finally re-write the 
equations in the following unique form: 



J,XX 



, M 2 ^(lj + 1) 



*i = 0, (54) 



where j — 0, 1, 2., and 



4 = £, £ = l,2,... 
£ x = £ + 1, £ = 0,1,2,. 
£ 2 = £-1, £=1,2,... 



(55) 
(56) 
(57) 



Repeating analysis of the section II, we can easily see 
that the intermediate late-time behavior will be 



* ~i~ (£+3/2) sinr77i, £=1,2,... (58) 
~f-^+ 5 / 2 ) s inmf £ = 0,1,2,... (59) 
§ 2 ^ t-( £ + 1 /2) sinmt £ = 1,2,... (60) 

B. Asymptotic late-time behavior 

If in the regime x 3> 1, we take into consideration sub- 
dominant terms of order u/x in Eq. (|42|) - (|44p . we come 
to to the following equations: 
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1W + 



w z 1 + 



2 / J 

1 + - 



A 2 " 1 



*o + ^*i =0 



(61) 



*2.^ + 



w 2 IH 




-t(- 


■9 


A 2 










a; 2 



*2 + ^5-*l =0 



(62) 



w 2 fn 






■=) 


A 2 + 2" 










x 2 j 



*1 



V2A 



(*o + #2) = 0. 



(63) 



Following [28], one can introduce a new variable z = 
2\fm 2 — La 2 r, then one can see that the system of equa- 
tions (|6U63[) has a solution [28]: 

*j = C^ M (4 3 =0,1, 2. (64) 
where W^, is the Whittaker function [2^] 

q = (2lu 2 - m 2 )M/\/m? -co 2 , (65) 

and for p, as was shown in [28], one has the following 
system of equations 

(p 2 - A 2 - (1/4))C + v^ACi = 0, (66) 

(p 2 - A 2 - (1/4))C 2 + V2ACi = 0, (67) 

(p 2 - A 2 - (9/4))Ci + \/2A(C 2 + C ) = 0. (68) 

Therefore, one immediately has values of p: 

p = I + (1/2) + a, a = 0,±l. (69) 

As was shown in [17| for massive scalar field, the 
asymptotically late-time behavior of the perturbations 
governed by the wave equation which looks at asymptot- 
ically large r as a Whittaker equation, the decay law is 
universal and does not depend on I. Thus, in a similar 
fashion with [13], we find that the decay at asymptoti- 
cally late times t — > 00 again is 

*j ;~ r 5/6 sinmt, j = 0,1,2. (70) 

IV. CONCLUSIONS 

In the present paper we have considered evolution of 
perturbations of massive vector field at late times. For 
spherically symmetric case, the perturbation equations 
can be reduced to a single wave-like equation, and there- 
fore can be analyzed numerically with characteristic in- 
tegration method. In this case, time-domain picture of 
a signal has good agreement with analytical predictions 
for the considered range of times. 



General vector field perturbations can be represented 
by three scalar functions ^>i, (i = 0,1,2). At interme- 
diate late times r « ! « M/(Mm) 2 , these three func- 
tions decay according to different laws which depend on 
a multipole number £, * ~ t~( £ + 3 / 2 ), vf^ ~ t^ ( - e+5/2 \ 
^2 ~ t~( i+1 / 2 \ On the contrary, at asymptotically 
late times t — > 00, the decay law is universal, i.e. the 
same for all three functions and does not depend on I: 
\Fj ~ £~ 5 / 6 sin rat. Note again, that the same asymptot- 
ically late-time decay low was previously found in [2l| . 
[l7T ] for scalar and spinor (Dirac) perturbations around 
Schwarzschild black hole. Summarizing this we can con- 
clude that at asymptotically late times t — > 00, the decay 
low is universal and does not depend on spin of a field s 
and multipole number I: 



* s - t~ 5/(i smmt, s = 0,1/2,1. (71) 



The latter law is obtained analytically and confirmed 
numerically with high accuracy. 

For Schwarzschild-de Sitter black hole it is observed 
quite non-trivial tail behavior. In the considered asymp- 
totically de Sitter geometries, the power-law oscillatory 
tails are replaced by purely exponential tails, if m is 
small, or by oscillatory quasinormal decay, if the parame- 
ter m is large enough. Although there is some universal- 
ity for the tail behavior in the context of Schwarzschild 
geometry, this is not so in non-asymptotically flat back- 
grounds. 

An interesting question, which was beyond our study, 
is whether the asymptotically flat background universal- 
ity of tail behavior keeps in higher dimensions, and for a 
more general then Schwarzschild backgrounds? Yet, for 
Kerr black hole there is no hope to separate variables 
in the perturbed Proca equations, while for a Reissner- 
Nordstrom black hole such a separation is certainly possi- 
ble for spherically symmetric perturbations, and charged 
Proca field should be considered instead. We hope that 
future investigations will clarify all these points. 
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